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^SJ , Abstract. We present new interferometric observations of the continuum emission at mm wavelengths in the 

00 ' Seyfert galaxies NGC 1068 and NGC 3147. Three mm continuum peaks are detected in NGC 1068, one centered 

0^ ' on the core, one associated with the jet and the third one with the counter-jet. This is the first significant detection 

, of the radio jet and counter-jet at mm wavelengths in NGC 1068. While the fluxes of the jet components agree 

with a steep spectral index extrapolated from cm- wavelengths, the core fluxes indicate a turnover of the inverted 
cm- into a steep mm-spectrum at roughly ^50 GHz which is most likely caused by electron-scattered synchrotron 
. emission. As in NGC 1068, the spectrum of the pointlike continuum source in NGC 3147 also shows a turnover 

CL(' between cm and mm-wavelengths at ~25 GHz resulting from synchrotron self-absorption different to NGC 1068. 

JL I This strongly resembles the spectrum of Sgr A*, the weakly active nucleus of our own galaxy, and M81*, a link 

■ between Sgr A* and Seyfert galaxies in terms of activity sequence, which may display a similar turnover. 

-4— » 

, Key words, galaxies: active - Seyferts: individual: NGC 1068 & NGC 3147 - Radio continuum: galaxies 

> 
• I— I 

^ 1. Introduction NGC 1068, at a redshift of z = 0.004, is one of the 

H best studied Seyfert galaxies. It is regarded to be the 

. P_ , NGC 1068 and NGC 3147 are both part of the NU(clei archetype for unification schemes of Seyferts. Khachikian 
of)-GA(laxy) project (Garcfa-BuriUo et al. 2003) that re- ^ Weedman (1974) first classified NGC 1068 as Seyfert 2 
cently succeeded m mapping the molecular gas emission ^j^^ detection of unpolarized (direct-path) emission 

of CO(l-O) and C0(2-l) m 12 nearby Seyfert and ^-^^^ ^^^^ ^-^^^^ originating in the Narrow-Linc- 

LINER galaxies with high angular resolution and high sen- ^^^^^^ (^LR) . The NLR appears to be extended to the 
sitivity using the IRAM Plateau de Bure Interferometer North-East with respect to the nucleus (e.g. Groves et al. 
(PdBI). NUGA primarily aims at analysing the dynamics 2004). The radio jet is found to be located in the same 
and distribution of the gas on the smallest possible scales direction (e.g. GaUimore et al. 2004). Emission from the 
to disentangle the accretion mechanisms of the gas onto ^^^^^d line region (BLR) could also be detected through 
the nucleus which are supposed to be strongly correlated polarized (scattered-path) emission lines with broad line 
with the observed nuclear activity. As a by-product of this ^-^^^^ (Antonucci & Miller et al. 1985). This strongly 
project, information on the continuum emission at 1 mm g^gggsts that NGC 1068 contains a hidden Seyfert 1 nu- 
and 3 mm was also gained for NGC 1068 and NGC 3147. ^j^^^ supporting the unified theory for Seyfert galaxies. 
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The position of the nucleus could be precisely derived by 
radio continuum observations carried out by GaUimore et 



e-mail: krips@phl.uni-koeln.de 

* rr^v- u J u 4-- -4-1, 4-1, TDAA/r al- ( 2004 ) . Thcy Tcport a compact componcut With a flat 

ihis paper was based on observations with the IRAM \ / j t- t- t- 

Interferometer. IRAM is supported by INSU/CNRS (France), spectrum lying in a maser disk that is most likely the 

MPG (Germany) and IGN (Spain). central engine. NGC 1068 is one of the few (i.e. ~ 30) 
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Fig. 1. Continuum maps of NGC 1068 at 3mm {left panel) and 1 mm {right panel). Contour levels are from 3(7=1.1 to 
25.3 mjy/beam in steps of 3f7 at 3 mm and from 3(7=2.1 mJy/beam to 14.7 mJy/beam in steps of 2fT at 1 mm. The 
beam size at both wavelengths is indicated as filled black ellipse in the lower left corner of the maps: 3mm - 2.1" x 1.2", 
PA=32°; 1mm - 1.0"x0.6", PA=36°. The (0,0) position of the map corresponds to the radio core position (see caption 
of Table HJ. 
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Fig. 2. Continuum maps of NGC 3147 at 3mm {left panel) and 1 mm {right panel). Contour levels are from 3(7=0.72 to 
4.8 mJy/beam in steps of 3(7 at 3 mm and from 3(7=1.4 to 3.5 mJy/beam in steps of Icr at 1mm. The beam size at both 
wavelengths is indicated as filled black ellipse in the lower left corner of the maps: 3mm - 1.4"xl.2", PA=57°; 1mm - 
0.7" X 0.5", PA=60°. The (0,0) position of the map corresponds to the radio core position (see caption of Table P). 



extragalactic systems in which maser emission could be detected. With the aid of H2O maser emission, Greenhill 
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Table 1. Continuum parameters for the 3 mm and 1 mm continuum in NGC 1068 and NGC 3147. Parameters were 
determined by fitting Gaussian profiles to the respective components. 



NGC 1068 


Component 


RA 


Doc 


flux 


Size 




(offset)" 


(offset)^ 


density'* 


maj.xmin., PA 




["] 


["] 


[mjy] 


["x'M 


core-3mm 


O.liO.l 


0.2±0.1 


36±0.4 


(0.8±0.1) X (0.6±0.1), 70±10 


jet-3mm 


2.5±0.3 


4.1±0.3 


17±0.5 


(2.6±0.1) X (1.2±0.1), 39±2 


counterjet-3mm'' 


-1.4±0.5 


-3.8±0.5 


7±0.6 




core- 1mm 


0.1±0.2 


0.2±0.2 


22±0.8 


(0.6±0.1) X (0.47±0.04), 40±20 


jet- 1mm'' 


~3 


~5 


6±1.0 




NGC 3147 


core-3mm 


O.OiO.l 


O.liO.l 


5.3±0.2 


O.liO.S'^ 


core-lmm 


0.1±0.2 


-0.2±0.2 


2.8±0.5 


0.1±0.3= 



offsets with respect to qj2ooo=02''42'"40.71= and (5j2ooo=-00°00'47.94" for NGC 1068 (corresponds to the position of the 
radio core; Muxlow et al. 1996) and to aj2ooo=10''16'"53.65'' and <5j2ooo=73°24'02.70" for NGC 3147 (e.g. Ulvestad & Ho 
2001; Krips et al., in prep.). 
^ The SNR is not sufflcicnt to estimate a size, 
emission most likely point-like 

values have been determined by fitting Gaussian profiles to the visibilities. The uncertainties do not include systematic 
errors from the fiux calibration (~10-20%). 



& Gwinn (1997) determined the mass of the AGN to 
1.7x10^ MqIu the central 2.2 pc. The bolometric lumi- 
nosity of NGC 1068 can be derived to a few ~10~^Ledd 
indicating - in terms of Eddington luminosity - a slightly 
higher radiation efficiency of the accretion and/or a higher 
accretion rate than in other LLAGN. The molecular gas 
in NGC 1068 was mapped by Schinncrcr ct al. (2000, and 
references therein). It is distributed in a two arm spiral, 
along a bar and in a warped nuclear ring with two emis- 
sion peaks which arc interpreted as orbit crowding points. 
Schinnerer et al. (2000) report also on a detection of the 
3 mm continuum emission in NGC 1068: two components 
become visible at this wavelength, a 16(t peak at the radio 
core position and a 3cr peak at the position of the radio 
jet. Going to longer wavelengths, more continuum struc- 
ture appears (e.g., Ulvestad ct al. 1987, Gallimorc et al. 
1996 & 2004): the central continuum peak splits up into 
four components, of which two show steep radio spectra 
(i.e., /,y oc 1/" with a « —1) while the remaining two have 
flat to inverted spectra (i.e., a ~ —0.2 and a ~ 0.3). The 
inverted spectrum feature (labeled SI in Gallimore et al. 
1996) is associated with the nucleus of NGC 1068 whereas 
the remaining three components are supposed to be part 
of the extended jet and counter-jet. 

NGC 3147, at a redshift of z 0.009, is classified as 
Sbc galaxy hosting a Seyfert 2 nucleus (Ho et al. 1997). 
ROSAT observations (Roberts & Warwick 2000) unveil an 
X-ray source at the position of the nucleus which is associ- 
ated with a supermassive black hole. Assuming the empir- 
ical correlation between black hole mass and bulge stellar 
velocity dispersion (e.g., Gebhardt et al. 2000, Ferrarese 
& Merritt 2000), the black hole mass can be estimated to 
~4-10^ M0(Ulvestad & Ho 2001). The bolometric lumi- 



nosity can be determined to ^10~'*Lcdd (Ulvestad & Ho 
2001) indicating a slightly higher radiation efficiency of 
the accretion than in other LLAGN similar to NGC 1068. 
VLBA and VLBI radio observations also show a pointlikc, 
non-thermal continuum source at the position of the nu- 
cleus (Ulvestad & Ho 2001, Anderson et al. 2004, Krips et 
al. in prep.). The radio spectrum is inverted and no signs 
of variability are found on timescales of a few years. The 
CO emission is astonishingly similar to NGC 1068. Besides 
a double-peaked nuclear ring observed in C0(2 1), spiral 
arms and/or a ring- like structure become visible at larger 
scales in IRAM PdBI observations (Garcfa-Burillo et al. 
2005, Leon et al., in prep.). 

2. Observations 

2.1. NGC 1068 

Observations of NGC 1068 at 3 mm and 1 mm were car- 
ried out with the IRAM PdBI in February 2003 using all 
six antennas in A configuration. The bandpass was cali- 
brated on NRAO150 while phase and amplitude calibra- 
tion were performed on 0235-^164 and 0238-084. CRL618 
and MWC349 were chosen to determine the absolute flux 
scale. The receivers were tuned to the redshifted ^'^C0(1- 
0) line (i.e., 109.8 GHz) at 3 mm and to the redshifted 
^^00(2-1) line (i.e., 231.1 GHz) at 1 mm. A total band- 
width of 580 MHz with a spectral resolution of 1.25 MHz 
was used per sideband. Applying uniform weighting^ in 

^ Actually, when speaking of uniform weighting, it rather cor- 
responds to robust weighting which is a compromise between 
natural and pure uniform weighting. This is a convention in 
the IRAM GILDAS Packages. 
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Fig. 3. SED of the continuum emission in NGC 1068. 
Upper panel: 230 GHz, 115 GHz and 110 GHz fluxes of 
the core [filled circles), the jet (open black circles) and 
the counter-jet {open grey circles) compared to 22 GHz 
and 5 GHz VLA observations {triangles; GaUimore et al. 
1996). The filled triangles correspond to fiuxes still includ- 
ing contributions from the core-jet components visible on 
VLA scales (< 1"; Gallimore et al. 1996). Middle panel: 
Core-jet components derived by Gallimore et al. (1996) 
and extrapolated to mm wavelengths. Lower panel: Pure 
core spectrum. The filled triangles show the cm flux ob- 
tained by Galhmore et al. (1996; their SI component). 
The filled circles represent the (estimated) pure core flux 
at mm wavelengths for which the core-jet flux was sub- 
tracted (see text for more details). The filled grey circles 
indicate the uncorrected mm-fluxes. The dashed-dotted 
curves (ESSE) show the uncertainties of the fit. The thick 
grey line is the thermal FFA model. 



I PdBI (this paper) 

I MERLIN (Krips et ol., in prep.) 

1 VEBA (Ulvestad & Ho 2001; Anderson et al. 
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Fig. 4. Spectral energy distribution of the core emission in 
NGC 3147. The 230 GHz and 110 GHz fluxes of the core 
(filled circles; this paper) are compared with the VLBA 
2.3 GHz & 8.4 GHz fiux (Ulvestad & Ho 2001), the VLBA 
15.4GHz, 22.2 GHz & 43.2 GHz flux (Anderson et al. 2004) 
and the MERLIN 5 GHz and 1.8 GHz flux (Krips et al., in 
prep.). The dotted-dashed curves indicate the uncertain- 
ties of the fit. 



the mapping process, beamsizes are derived to 2.1x1.2" 
with a position angle (PA) of 32°at 3mm and to 1.0x0.6" 
with PA=36°at 1mm. 



2.2. NGC 3147 

i2CO(l-0) and 12C0(2-1) emission from NGC 3147 were 
observed with the IRAM PdBI in February 2002 and 2004 
using all six antennas in A, B, C, and D configuration. 
3C345 was used as bandpass calibrator, while 1546-1-027, 
1044-1-719 and 0838-1-710 were chosen as amplitude and 
phase calibrators. The flux was calibrated on MWC349 
and 3C273. The receivers were tuned on the redshifted 
i2CO(l-0) line (i.e., 114.2 GHz) at 3 mm and to the red- 
shifted i2C0(2~l) line (i.e., 229.9 GHz) at 1 mm. The to- 
tal bandwidth was set to 580 MHz in each sideband with 
a spectral resolution of 1.25 MHz. The natural beamsizes 
are determined to 1.4" x 1.2" with PA=57°at 3 mm and to 
0.7" X 0.5" with PA=60°at 1 mm. 



3. Results 

The continuum maps at 3 mm and 1 mm were computed 
from those channels in which no line emission is expected 
(i.e.. Upper Side Band (USB) at 3 mm and Lower Side 
Band (LSB) at 1 mm) . Continuum emission was detected 
at 3 mm and at 1 mm in both sources. 

All cm and mm fluxes are accurate at least to within 
~ 10 — 20% accounting also for uncertainties in the flux 
calibration and systematic errors. 
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Fig. 5. Spectral Energy Distribution of NGC 1068 and 
NGC 3147 from radio wavelengths to X-rays. Radio and 
mm fluxes are from this paper while the remaining values 
are taken from NED and Ulvestad & Ho 2001. 

3.1. NGC 1068 
3.1.1. 3 mm 

The 3 mm continuum map is shown in Fig. ^ revealing 
three different components, one associated with the nu- 
cleus, one with the jet and one with the counter-jet. The 
jet component is elongated in North-East direction, agree- 
ing with results found at cm wavelengths. Also a 6a fea- 
ture is detected on the opposite side of the nucleus indi- 
cating the counter-jet. Schinnerer et al. (2000) detected 
the core with a flux of 36±5 mJy and the jet with 3cr to 
15±5 mJy at a slightly higher frequency (115 GHz) but 
missed the counter jet. This agrees well in flux and posi- 
tion with our new data (see Table 

Three elliptical Gaussian profiles were fitted to our 
3 mm data. The results are listed in Table Q The po- 
sition of the core is determined at aj20oo= 02M2™40.71" 
and (5j2ooo=~00°00'47.7" consistent with the core posi- 
tion derived at 6 cm (Gallimore et al. 2004). The jet and 
counter-jet fluxes are plotted in Fig. |2| at different fre- 
quencies including radio fluxes estimated by Gallimore et 
al. (1996). The diagram shows a steep spectral index of 
a « —0.5 {fi, oc ly") for the jet and for the counter-jet. 
This is consistent with optically thin synchrotron radia- 
tion that is expected for jet emission. The situation for the 
core is more complex since VLA maps (Gallimore et al. 
1996) unveil further subarcsecond components which can- 
not be resolved with the IRAM PdBI: The PdBI "core" 
is resolved by VLA into two further jet features (NE,C) 
with steep spectra and two flat to inverted components 
(S1,S2) one of which is associated with the nucleus (la- 
beled SI in Gallimore et al. 1996). An estimate from 
lower resolution radio maps (Gallimore et al. 1996), not 
distinguishing between the respective subarcsecond com- 
ponents, indicates a steep spectrum for the overall (in- 
ner '^2") region we designated as "core(-l-corejet)" (i.e. 
S1-I-S2+NE4-C; a=-0.9; upper panel of Fig. ^. Thus, 
optically thin synchrotron radiation from the subarcsec- 



ond jet (=core-jet; NE-I-C-I-S2) most likely dominates the 
emission at mm- wavelengths. However, since the fluxes of 
the subarcsecond radio (22 GHz) spots (namely NE,C and 
S2 in Gallimore et al. 1996) are known, their contribution 
to our "core(-l-corejet)" mm-component can be estimated 
to separate it from the nuclear component SI and at least 
a rough estimate for SI at 3 mm can be given. Assuming 
the spectrum calculated from Gallimore et al. (1996; mid- 
dle panel in Fig. ^ for these features, we find that ^--^60% 
of the 3 mm core(-l-corejet) still comes from these com- 
ponents (middle panel of Fig. |3J|. Subtracting this value 
from the 3 mm fluxes (115 and 110 GHz) allows to com- 
pare both with the core (SI) fluxes at 5 GHz of 12±1 mJy 
and at 22 GHz of 19.1±0.6 mJy (VLA, Gallimore et al. 
1996). The lower panel of Fig.^lindicates a drop of the fiux 
from the power law spectrum with a = 0.3 towards the 
3 mm data. Accounting for the narrower beam at the cm- 
wavelengths compared to the mm-observations and thus 
resolution effects, the drop might be even more extreme^. 
Another caveat to keep in mind is possible variability. 
However, the mm fiuxes give no hints for any variabil- 
ity between different epochs although a fiux variation of 
~30% has been reported for one of the core-jet compo- 
nents (C) (Gallimore et al. 2004). This needs to be verified 
with future observations. 

3.1.2. 1 mm 

In our PdBI data from 2003, we clearly detect 
(>20(T=22 mJy/beam) continuum emission at 1 mm 
(Fig. nj, while Schinnerer et al. (2000) report only an 
upper limit of the 1 mm continuum of 6 mJy. However, 
a refined analysis of the earlier Schinnerer et al. (2000) 
maps shows that also their data contain continuum emis- 
sion which is consistent to 10% with what we calculate 
from the new observations. The position of the core de- 
rived from the 1 mm data is consistent with the position 
derived at 3 mm. The fiux of the 1 mm core component fol- 
lows the power law fitted to the lower frequencies in Fig. El 
(upper panel). One has to keep in mind, however, that the 
core emission at 1 mm might also still be biased by the dif- 
ferent core-jet components (NE,C and S2; Gallimore et al. 
1996) within the central 1" detected with VLA and VLBA. 
Assuming the spectral indices estimated from the radio 
fluxes and taking the (VLA) fiuxes at 22 GHz (Gallimore 
et al. 1996), the contribution of these further components 
(NE,C, S2) to the pure core (SI) emission at 1 mm can 
be estimated to ~60% (middle panel of Fig. 0). The fiux 
of the estimated mm core-jet might be even higher ac- 
counting for the higher angular resolution of the cm data 
compared to the mm data and thus of resolution effects, 
so that the 1 mm flux of the pure core could still be lower. 
Subtracting the contribution of the core-jet from the 1mm 
flux and comparing it with the core flux derived at cm 
wavelengths by Gallimore et al. (1996), we find further 

^ Note that the beamsize at 5 GHz is however larger than at 
22 GHz. 
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support for a turnover between cm and mm wavelengths 
(lower panel of Fig. Even if we do not apply a core- 
jet correction and assume hence the 1 mm emission to be 
mainly dominated by the core alone, the 1 mm flux is al- 
ready significantly below the inverted power law spectrum 
produced by the cm and 3 mm fluxes. Although resolution 
effects between 3 mm and 1 mm have to be taken into ac- 
count, they should only play a significant role for the ex- 
tended core-jet components but not for the core itself since 
it appears to be quite compact at cm wavelengths (VLA 
scales). Thus, accounting for all arguments, we think that 
the turnover is a real intrinsic property of the NGC 1068 
nucleus and no artefact of resolution effects or calibration 
uncertainties. 

In general, three mechanisms can be proposed to 
explain such a low-frequency turnover in NGC 1068 
(e.g., Roy et al., 1998, Gallimore et al., 1997 & 2004): 
(1) synchrotron emission from a plasma with electron 
self-absorption (SSA), (2) synchrotron emission which is 
electron-scattered in the obscured central region (ESSE), 
or (3) thermal free-free absorption from an X-ray heated 
(thermal FFA), ionized plasma. The spectra in all three 
cases can be fitted with the following equation (e.g., Kraus 
1986): 



^ "thick 

U oc { —] I 1 - exp 



V 



where vq is in units of 10^ Hz, and athick and atinn are 
the respective spectral indices for the optically thick/thin 
part of the spectrum. In the first case of synchrotron self 
absorption (SSA) and also in the third case of thermal 
free- free absorption (FFA), /t^ is 1. In the second case 
of electron-scattered synchrotron emission (i.e. free-free 
absorption by a scattering ionised gas; ESSE), /t- is equal 
to exp(— T//(i/)) with: 



Tff{v) = ^ 0.08235 X r-i-35n2Ljy-2 i 



(2) 



Te is the electron temperature in units of Kelvin, rig is 
the mean square electron density in cm~^, and L is the 
path length along the sight line in pc (e.g., Altenhoff et 
al., 1960; Mezger & Henderson 1967). 

Gallimore et al. (1996) already excluded the first 
scenario based on the low brightness temperatures of 
~3xl0^ K derived at 5 GHz. They argue that this tem- 
perature is too low for synchrotron self absorption which 
generally becomes important at much higher brightness 
temperatures (> 10* K; e.g., Williams 1963, Kellerman 
& Owen 1989). We therefore consider only the two lat- 
ter cases in the following of either electron-scattered syn- 
chrotron emission (ESSE) or thermal free-free absorption 
(FFA) for the SI component in NGC 1068. As thermal 
free- free emission is found to have generally athick — 2 and 
CKthin = (e.g., Kraus 1986 (Orion or Rosette nebula); Roy 
et al., 1998), this does not appear to be consistent with the 
measured spectrum in which the data point at 1 mm lies 

^ /t is only an artificial parameter to summarize all three 
cases in one equation. 



significantly below the thermal FFA curve (see boldfaced 
grey solid curve in the lower panel of Fig.|2|). However, one 
important caveat to be considered are resolution effects 
which might result in an artificial drop of the spectrum 
between 3 mm and 1 mm. This cannot be totally excluded 
although it seems likely that this effect has a stronger in- 
fluence on the extended jet and core-jet components than 
on the rather compact core (SI). Assuming the ESSE case 
(see solid and dashed-dotted curves in the lower panel of 
Fig. 131) , the observed spectrum can be nicely reproduced. 
The fitted optical depth f is consistent with an absorb- 
ing plasma at Tg — lO^K, and an emission measure of 
r?^L — 10^cm~^pc agreeing with the values used by Roy 
et al. (1998) and Galhmore et al. (2004). The turnover fre- 
quency can be derived to 50±10 GHz for NGC 1068. The 
error was calculated from the fit assuming equal weighting 
factors for the data points. 

Besides the core emission, a further weak (~5cr) con- 
tinuum peak is detected at (roughly) the position of the 
large scale radio jet. The flux is compatible with the steep 
jet spectrum (/,y cx i/""-^) extrapolated from lower fre- 
quencies. 

3.2. NGC 3147 



/r,(l) 3.2.1.3 mm 



The 3 mm continuum appears pointlike consistent with ra- 
dio observations at different angular resolutions (VLBA: 
Ulvestad & Ho 2001, VLBI & MERLIN: Krips et al., in 
prep.). While the radio fluxes span an inverted spectrum 
with a w 0.1 — 0.4 (/^ cx i/"), the 3 mm continuum sig- 
nificantly falls below the spectrum extrapolated to the 
mm domain from radio frequencies (Fig. 2J. This can- 
not be explained by resolution effects. First, NGC 3147 
is pointlike almost on all angular scales which is sup- 
ported by identical radio fluxes at one frequency but de- 
rived from different angular resolutions. Second, the radio 
fluxes derived with MERLIN (i.e., beamsizes are <0.2") 
are higher than the mm-flux determined with the IRAM 
PdBI (beamsize is ^1"). Variability equally appears to be 
unlikely. Observations at different epochs, different wave- 
lengths (18cm, 6cm, 3mm & 1mm), and with different in- 
struments (IRAM PdBI, MERLIN & VLBI) give no hints 
for any variation above 10% of the flux in NGC 3147, nei- 
ther at cm nor at mm wavelengths (Ulvestad & Ho 2001, 
Krips et al., in prep.). 

3.2.2. 1 mm 

As already expected from 3 mm and cm wavelengths, the 
continuum emission at 1 mm is pointlike as well. Based 
on the results of the radio data, resolution effects can 
thus be almost discarded. The flux derived at 1 mm fur- 
ther decreases compared to 3 mm and higher wavelengths 
(see Fig. 01). Thus, the turnover in the spectrum from cm 
to mm wavelengths indeed is an intrinsic property of the 
active nucleus in NGC 3147. The relatively high bright- 
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ness temperature derived at 5 GHz of at least ~ 10^ K 
(e.g., Anderson et al. 2004, Krips et al., in prep.), sug- 
gests synchrotron self-absorption to be responsible for the 
turnover in NGC 3147, different to NGC 1068. By fitting 
Equation ^ to the data, the turnover frequency can be 
derived to 25±10 GHz. 

4. Correlation between radio/mm and X-rays 

Fig. El presents the overall spectrum of NGC 1068 and 
NGC 3147 from radio to X-ray frequencies. For complete- 
ness, infrared and optical fluxes are plotted additionally. 
The cm/mm data points correspond to those shown in 
Fig. (bottom) and Fig. H 

Through Synchrotron Self (inverse) Compton scatter- 
ing (SSC), the submm- wavelength photons can be up- 
scattered into the NIR and X-ray domain. Using the for- 
mula given by Marscher (1983) and Gould (1979), we can 
constrain the size of the emitting regions with the ob- 
served X-ray fluxes in NGC 1068 (in the case of ESSE) 
and NGC 3147. To reproduce the observed values, a radius 
of the order of ~100i?s {Rs is the Schwarzschild radius) 
for the emitting region in NGC 1068 and NGC 3147 are 
needed, assuming a Lorentz factor 7e = (1 — v/c)~^^''^ of a 
few thousand. While this is consistent with the radio data 
of NGC 3147, it appears to be too small for NGC 1068, rel- 
ative to the VLBI images of Gallimore et al. (2004), which 
suggest a size 10 times larger. Although the X-ray source 
is assumed to be of the same order, the radio emission 
originates from an area which is at least an order of mag- 
nitude larger (Gallimore et al. 2004). This indicates that, 
at least in the case of NGC 3147, the underlying physical 
mechanisms producing the radio and X-ray emission are 
correlated with each other through SSC whereas the sit- 
uation in NGC 1068 might be more complicated due to 
the scattering process and the equation used above might 
not be appropriate anymore. The correlation between X- 
rays and radio fluxes via SSC was already reported for 
quasars (e.g., Roberts & Warwick 2000), the equivalents 
to NGC 1068 and NGC 3147 at a much higher activity 
level, and was recently also found for low-power accreting 
black holes (e.g., Falcke et al. 2004). As a thermal FFA 
cannot be totally excluded for NGC 1068, a correlation 
between the radio continuum and X-ray emission is also 
expected in this case since the absorbing ionized plasma 
producing the thermal free-free emission is heated by X- 
rays coming from the AGN. 

The turnover frequency of a synchrotron spectrum de- 
pends on the magnetic field strength among others. Taking 
again the equation used by Gould (1979) and Marscher 
(1983), the magnetic field density can be estimated in the 
case of NGC 3147. However, this estimate also involves 
an assumption of the source size and is very sensitive to 
the assumed value so that it must be taken with caution 
and rather understood as a consistency check. However, 
assuming SSC mechanism as discussed above, the X-ray 
flux set a constraint on the size of the emitting region 
for NGC 3147 yielding a radius around 100i?s. If we as- 



sume this size for NGC 3147, then the turnover frequency 
of 25 GHz at a flux of 12 mJy implies a magnetic field 
strength of ~10 Gauss. This is ~3 times lower than the 
magnetic field of Sgr A*, for an accretion rate below the 
Bondi value (Quataert & Gruzinov 2000, Eckart et al. 
2004). 

5. Summary &. Discussion 

Continuum emission was detected in NGC 1068 and 
NGC 3147 with the IRAM PdBI at 3 mm and 1 mm. 
While the maps of NGC 1068 unveil several extended com- 
ponents in form of a jet, a counter-jet and an extended 
core, the continuum emission in NGC 3147 appears point- 
like. The jet and counter-jet in NGC 1068 were detected 
for the first time with a large signal-to-noise ratio at mm- 
wavelengths. Their positions and fluxes are consistent with 
values obtained from radio observations. Extrapolating 
the influence of the core-jet (<!") components seen in 
VLA maps (Gallimore et al. 1996), i.e. at higher angu- 
lar resolution, on the measured mm-fiux at the position 
of the core, a turnover in the pure core spectrum is in- 
dicated at roughly ^50 GHz. However, further observa- 
tions at higher angular resolution are needed to separate 
the different components within 1" at mm- wavelengths. 
This will allow to strengthen the case for the turnover in- 
dicated here. The situation in NGC 3147 is significantly 
less speculative. The pointlike structure of its continuum 
emission discards the influence of resolution effects on the 
derived fluxes induced by different angular resolutions. 
The spectrum clearly reveals a turnover from flat (cm) to 
steep (mm). The turnover frequencies can be determined 
to -50 GHz for NGC 1068 and to -25 GHz for NGC 3147. 
While the turnover in the spectrum of NGC 1068 is most 
likely caused by electron-scattered synchrotron emission 
(although thermal free-free absorption cannot be totally 
discarded with our new mm-data), the turnover in the 
spectrum of NGC 3147 can be associated with synchrotron 
self-absorption. 

The spectra of NGC 1068 and NGC 3147 both remind 
of the one of Sgr A* at the center of our own galaxy 
(e.g. Meha & Falcke 2001). Also in M81'^, the center 
of M81, hints have been found that its spectrum might 
have a turnover between cm and mm- wavelengths. The 
bolomctric luminosity (lO^^Lcdd, Renter & Lesch 1996) 
of M81* is intermediate between that of Sgr A* (— 10~^- 
lO^i^Lcdd, Meha k Falcke 2001; Yuan et al. 2003) and the 
2 Seyfert galaxies NGC 1068 (-lO-^Lodd) and NGC 3147 
(— 10~*Lodd) and can thus be regarded as link between our 
weak galactic center and low-luminosity AGN. The similar 
behaviour of their spectra might be a further hint for this. 
Turnovers in the synchrotron spectra of compact radio 
sources have been already identifled in quasars before, the 
other extreme in terms of activity (e.g., Lobanov 1998, 
Bloom et al. 1999). Giga-Hertz-Peaked (GPS) quasars 
with turnover frequencies at a few GHz are known since 
quite a while but quasar spectra peaking at several 10 GHz 
have been found as well (e.g.. Bloom et al., 1999). This fur- 
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ther connection from AGN with low radiation efficiencies 
of tlie accretion and/or higlily sub-Eddington accretion 
rates over those at somewhat higher efficiencies/rates to 
the most extreme cases with accretion rates close to the 
Eddington limit allows to set tighter constraints on the 
models suggested to explain the different activity levels. 
Sgr A* in particular provides the unprecedented opportu- 
nity to witness the accretion mechanisms at play on scales 
which are impossible to reach for extragalactic AGN but 
there is hardly any accretion going on in this case. Studies 
of the galactic center stimulated the re-emergence of radia- 
tively inefficient accretion flow models (RIAFs, see Mclia 
& Falcke 2001; Quataert 2003), including the advection 
dominated accretion flow approach (ADAF, e.g., Narayan 
& Yi 1994) and its derivatives. These models are char- 
acterized by low efficiency in converting thermal energy 
into electromagnetic radiation, but differ in terms of the 
inner boundary conditions such as the presence of strong 
winds (e.g., Blandford & Begelman 1999) and/or convec- 
tion (Quataert & Gruzinov 1999), each of which lead to 
differing predictions of the detailed spectral energy distri- 
bution (SED) and its variability. M81* will be an ideal 
testbed to broaden the available observational data base 
and verify, i.e. further constrain the existing models for 
sub-Eddington accretion. This will enable to set a closer 
link to somewhat higher or more efficiently radiating sys- 
tems like NGC 1068 and NGC 3147. Anderson et al. (2004) 
confronted the measured SED of NGC 3147 with a pre- 
dicted spectrum from ADAF and jet models. The simple 
ADAF models reproduce most of the observed spectrum 
in NGC 3147 but result in a too low radio emission. The 
same is probably true for NGC 1068 when comparing its 
spectrum with their models. Including a compact radio 
jet to the ADAF model (i.e., the optically thick part of 
the base of a (relativistic) jet (see also Athreya et al. 
1997)), the spectra can be better simulated. This sup- 
ports recent studies which have even found evidence for 
the jet dominating the nuclear radio emission in LLAGN 
over RIAFs (e.g., Nagar et al. 2001, Falcke et al. 2004) 
while omission from black holes accreting close to the 
Eddington rate arc rather dominated by their accretion 
disks. NGC 1068 clearly shows jet emission on different 
angular scales. Although thermal free-free emission in the 
core cannot bo totally discarded with our now mm-data, 
the mm/cm emission appears to be rather consistent with 
synchrotron emission from a compact component in the 
core (e.g., the base of the jet) which is electron-scattered 
by the obscuring torus. NGC 3147 has so far given no hints 
for any extended jet down to a scale of ~lpc but a com- 
pact jet cannot be excluded. Although the black hole in 
NGC 3147 is ten times more massive than in NGC 1068, 
its bolometric luminosity is an order of magnitude less 
than in NGC 1068 indicating a lower radiation efficiency 
of the accretion. 
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